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A fluid jet with a finite angular velocity is subject to centripetal forces in addition to surface 
tension forces. At fixed angular momentum, centripetal forces become large when the radius of the 
jet goes to zero. We study the possible importance of this observation for the pinching of a jet within 
a slender jet model. A linear stability analysis shows the model to break down at low viscosities. 
Numerical simulations indicate that angular momentum is expelled from the pinch region so fast 
that it becomes asymptotically irrelevant in the limit of the neck radius going to zero. 



I. INTRODUCTION 



A fluid jet emanating from a nozzle will become unstable and break up due to surface tension. Some 30 years ago, 
a series of papers |]-|| investigated the modifications Rayleigh's classical analysis would undergo if the jet performed 
a solid-body rotation. Such a rotation is easily imparted by spinning the nozzle at the appropriate frequency. The 
somewhat surprising result of the linear analysis is that the rotation always destabilizes the jet, a wavenumber k being 
unstable if 
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^ 0<fcr o <(l + i- 1 ) 1 / 2 , (1) 
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. with 

L = 1 /{ptfrl). (2) 

O i' Here 7 is the surface tension, p the density, fl the angular frequency, and ro the unperturbed jet radius. Note that fl 
appears in the denominator, so no rotation corresponds to L — -> 00, for which the stability boundary < kr$ < 1 found 
by Plateau is recovered. The theoretical growth rates were found to be in reasonable agreement with experiment Q] 
. and growth of disturbances for Zero larger than 1 was confirmed. 

' Recently it was pointed out by Nagel Q that rotation might have an even more dramatic effect for the highly 
nonlinear motion near the point where the neck radius goes to zero. Assume for the sake of the argument that a 
cylinder of fluid of length w pinches uniformly, i.e. it retains its cylindrical shape. Then the total angular momentum 
is 



on : M = (3) 

and the volume V = ttt^w is constrained to remain constant as ro goes to zero. The total interface pressure corre- 
ct sponding to the outward centripetal force is found to be p c — pr ,2 ] fl 2 /2, and thus 

p c = 2M 2 /(V 2 pr 2 ). (4) 

^ , As ro goes to zero, this outward pressure will dominate the surface tension pressure 7/7*0, raising the possibility that 
rotation is a singular perturbation for pinching: an arbitrarily small amount of angular momentum, caused by a 
symmetry breaking, could modify the breaking of a jet. 

However, a jet does not pinch uniformly, but rather in a highly localized fashion ||. If the above argument is applied 
to pinching, it must correspond to a rapidly spinning thread of fluid surrounded by almost stationary fluid. Frictional 
forces represented by the viscosity of the fluid will lead to a diffusive transport of angular momentum out of the pinch 
region, thus reducing its effect. Determining which effect dominates requires a fully nonlinear calculation, including 
effects of surface tension, viscosity, inertia, and centripetal forces. In the spirit of earlier work for fl — we 
will derive a one dimensional model, which only takes into account the leading order dependence of the velocity field 
on the radial variable. This will be done in the next section, together with a comparison of the linear growth rates 
between the model and the full Navier-Stokes calculation. In the third section we analyze the nonlinear behavior. 
First we investigate possible scaling solutions of the model equations, then we compare with numerical simulations. 
In the final section, we present some tentative conclusions. 
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II. THE MODEL 



In our derivation of the slender jet model we closely follow M. The Navier-Stokes equation for an incompressible 
fluid of viscosity v read in cylindrical coordinates: 

dtv r + v r d r v r + v z d z v r — v^/r = —d r p/p + v(d 2 v r + d 2 v r + d r v r /r — v r /r 2 ), 
d t v z + v r d r v z + v z d z v z = -d z p/p + v(d 2 r v z + d 2 v z + d r v z /r), 

d t v^ + VrdrVcf, + v z d z v^ + v r v$/r = v(d 2 v$ + d z Vj, + d r v^,/r - v^/r 2 ), 
with the incompressibility condition 

d r v r + d z v z + v r /r = 0. 

Here we have assumed that the velocity field does not depend on the angle </>. Exploiting incompressibility, v z and v r 
can be expanded in a power series in r: 

v z (z, r) — v (z) + V2(z)r 2 + . . . (5) 
vr(z,r) = —^r ~ -^r 3 - . . . . 

Here and in the following a prime refers to differentiation with respect to z. 

The crucial trick to make an expansion in r work in the presence of rotation is to rewrite v$(z, r) in terms of the 
angular momentum per unit length £(z) of the corresponding solid body rotation: 

v*{z,r)= 7T 2 p ^l ) r + br3 + ---- (6) 

Here h(z) is the local thread radius, hence no overturning of the profile is allowed. Just as without rotation, the 
equation of motion for h(z,t) follows from mass conservation based on the leading order expression for v z : 

d t h + v a h' = -v' h/2 . (7) 

Finally, the pressure is expanded according to 

p(z,r) = po(z) + p 2 (z)r 2 + . . . . (8) 

Plugging this into the equation of motion for v r , to leading order in r one finds the balance 

2i 2 

P2 = «w>' (9) 

while the leading balance for the w z -equation remains 

d t v + v v' = -p' /p + v{Av 2 + <). (10) 

Lastly, the i^-equation leads to 

d t £ + £v' + Uvoti/h + v a h 4 (e/h 4 )' = vh A {^pb + (£/h 4 )") (11) 

to leading order. 

Equations (|9|)-([ll]) contain the unknown functions po, v 2 , and b which need to be determined from the boundary 
conditions. The normal stress balance nun = 7K gives 

Po +P2h 2 =jk- v' , 

where n is the sum of the principal curvatures. As in the case without rotation, the tangential stress balance next = 
gives 
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-3v' Q h' - v'^h/2 + 2v 2 h = 
for t pointing in the axial direction, but a new condition 

nphb = h'(£/h 4 )' 

for t pointing in the azimuthal direction. Putting this together, one is left with a closed equation for h, Vq, and £: 

dth + vti = -v'h/2 (12) 
d t v + vv' = + -^-(£ 2 /h 6 )' + 3v(v'h 2 )'/h 2 

d t i + (viy = v{h i {i/h i ) i y, 

where we have dropped the subscript on vq. The same equations were derived independently by Horvath and Hubcr 
& 

The most obvious way to test this model is to compare with the known results for the stability of the full Navier- 
Stokes equation. To that end we linearize (Oh about a solution with constant radius ro and rotation rate fi: 



h(z,t) = r Q (l + ee wt cos(fcz)) (13) 

v(z,t) = -2e~e w * sinfe) 
k 



t{z,t) = |pOr^(l + ee"'acos(fcz)). 



Eliminating a, this leads to the equation 
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3 4fc 2 2 k 2 - 2 x R jT2 /D 2-- 



— w 2 + — (-1 + P + 7T 1 + &k 2 /Re 2 )uj + -— (-1 + P - 7T 1 ) = 0, (14) 
ite 2 2ite 



where k = kro and uj — uj(prQ) 1 ^ 2 /j 1 ^ 2 are dimensionless. We have introduced the Reynolds number Re = 
(t^o) 1 ^ 2 /(p 1 ^ 2 ^), based on a balance of capillary and viscous forces. Note that this convention differs from that 
of 0]. Putting uj = one reproduces the exact stability boundaries (|l|). However one can see that the inviscid limit 
Re — > oo is a very singular one, in disagreement with the full solution. Namely, for this limit one finds the three 
branches 

^fa-p-i- 1 ), ^=^ p g:f + + /_; 1} - as) 

Thus 0)3 is the only unstable mode in the range 1 — L^ 1 < k 2 < 1 + but goes to zero when the viscosity becomes 
small. The reason for this behavior, which is not found in the solution of the full equations, lies in the appearance of 
a very thin boundary layer for small viscosities Namely, Rayleigh's stability criterion for a rotating fluid implies 
that the interior of the fluid is stabilized. This forces any disturbance to be confined to a boundary layer of thickness 



2Qk 



near the surface of the jet, and S becomes very small for k « 1. But this additional length scale is not captured by 
our slender jet expansion. Only for high viscosities is the boundary layer smoothed out sufficiently, and from (|lj) one 
finds the dispersion relation 

ca = ^(l-P + L- 1 ), (16) 
6 

which is consistent with the full theory in the limit of small k. 

III. NONLINEAR BEHAVIOR 

Our main interest lies of course in the behavior close to pinch-off. Close to the singularity, one expects the motion 
to become independent of initial conditions, so it is natural to write the equations of motion in units of the material 
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parameters of the fluid alone. In addition to the known units of length and time, £ v — v 2 p/"f and t v — v 3 p 2 /j 2 , 
one finds an angular momentum scale £q = v 5 p 2 /^ 2 . Note that this scale is only about 1.9 • 10 -14 g cm/s for water, 
corresponding to a frequency of 2 • 10 _11 s _1 for a 1 mm jet, so even the smallest amount of rotation will be potentially 
relevant. Rewriting the equations of motion ( |l2j ) in units of ^„,i„, and £q, one can effectively put p = v = 7 = 1, 
leading to a universal form of the equations, independent of the type of fluid. 
In addition, one can look for similarity solutions [[llj of the form 

h(z 1 t) = t' ai 0(z'/t'P) (17) 
v{z,t) = t' a2 ip(z'/t^) 

e(z,t) = t' a °x(z'/t' ls ), 

where t' — to — t and z' = z — zq are the temporal and spatial distances from the singularity where h goes to zero. 
We have assumed that everything has been written in units of the natural scales £„, t„, and £ . By plugging (|l7]) into 
the equations of motion, and looking for solutions that balance the i'-dependence, one finds a unique set of values for 
the exponents: 

a x = 1, a 2 = -1/2, a s = 5/2, (3 = 1/2. 

In addition, one obtains a set of three ordinary differential equations for the similarity functions <j), and x- So far 
we have not been able to find consistent solutions to these equations, which match onto a solution which is static on 
the time scale t' of the singular motion. This is a necessary requirement since the fluid will not be able to follow the 
singular motion as one moves away from the pinch point. 

This negative result is consistent with simulations of the equations for a variety of initial conditions. To avoid 



spurious boundary effects, we considered a solution of (12) with periodic boundary conditions in the interval [—1,1] 



and an additional symmetry around the origin. This ensures that the total angular momentum is conserved exactly. 
We took the fluid to be initially at rest and the surface to be 

/ii„if(*) = r (l + 0.3cos(27rar)), (18) 

with tq — 0.1. The angular momentum was distributed uniformly with the initial value £mit, corresponding to 
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n 2 ipr% 
4 £ 2 ' 
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Figures [l] and || show a numerical simulation of (|T^) with Re — 4.5 and L = 0.25 using a numerical code very 
similar to the one described in jL2|. Written in units of the intrinsic angular momentum scale, £i n u/£o — 6 • 10 3 , so 
£ is potentially relevant. The thread pinches on cither side of the minimum, pushing fluid into the center. As seen 
in the profiles of £, the angular momentum is expelled very rapidly from the regions where h goes to zero and also 
concentrates in the center. This is confirmed by a plot of the minimum of £ versus the minimum of h. On the basis of 
the similarity solution (|l7|), a power law £ m i n oc is to be expected. Instead, Fig. |^ shows that £ m in decays more 
rapidly, making angular momentum asymptotically irrelevant. Indeed, a comparison of the similarity function <fi as 
found from the present simulation shows perfect agreement with the scaling theory in the absence of rotation p"l| ] . The 
behavior of £ m in should in principal be derivable from the linear equation for £ with known time-dependent functions 
h(z,t) and v(z,t). Unfortunately, £ m in does not seem to follow a simple scaling law except below h = 3 • 10 -4 , where 
the power is close to 3.13. It is not clear how to extract this power analytically. 

One might think that by increasing the angular momentum the system would cross over to a different behavior. To 
test this, the initial angular momentum was doubled to give L — 0.0625. At L — 0.5 centripetal and surface tension 
forces are balanced, so decreasing L significantly below this value will cause rotation to be important initially. Indeed, 
instead of pinching down immediately, the fluid is first pulled into a narrow disc, while the radius of the surrounding 
fluid remains constant, cf. Fig. |j. Eventually this outward motion stops, as surface tension and centripetal forces 
reach an equilibrium. Only then does the fluid pinch down at the edge of the disc. The behavior close to the pinch 
point is however exactly the same as for smaller angular momentum. As a word of caution, one must add that our 
model is certainly not valid at the edges of the disk, where slopes become very large. In fact, the very sharp gradients 
encountered in this region may be due to the fact that the fluid really wants to develop plumes. As is observed for low 



viscosity 12 , the viscous terms prevent the interface from overturning, but at the cost of producing unrealistically 



sharp gradients. 
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FIG. 1. The height profile in a numerical simulation with Re = 4.5 and L = 0.25. Shown are the initial condition, and the 
times when the minimum has reached h m in = 10 -1 ' 5 , 10 -2 , and 10 -5 , at the end of the computation. 
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FIG. 2. The angular momentum profiles £(z,t) corresponding to Fig. |l]. 
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FIG. 3. The minimum value of the angular momentum as function of the minimum height. It is found that l m i n decreases 

5/2 

faster than h^^, which would exactly balance surface tension and centripetal forces. 
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FIG. 4. A numerical simulation with twice the angular momentum of Fig. |l|. The height profiles are shown in time intervals 
of 0.05 and at the end of the simulation. Centripetal forces draw the fluid out into a disc. 
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IV. CONCLUSIONS 



The present investigation is only a first step towards the understanding of the role of rotation in droplet pinching. 
A major challenge lies in finding a description valid at low viscosities. This can perhaps be done by incorporating the 
boundary layer structure near the surface into the slender jet approximation. The relevance of this lies in the fact 
that angular momentum is potentially more important at low viscosities, when there is less frictional transport out 
of the pinch region. In fact it can be shown that the inviscid version of (|l2|) does not describe breakup at all, since 
centripetal forces will always dominate asymptotically. This result is of course only of limited use since the model 
equations are definitely flawed in that regime. 

In addition, there remains the possibility that a region in parameter space exists where angular momentum modifies 
breakup even at finite viscosities. We cannot make a definite statement since the additional variable makes it hard 
to scan parameter space completely. Finally, spinning jets have not received much attention in terms of experiments 
probing the non-linear regime. The discs found at high spinning rates (cf. Fig. ^ are a tantalizing new feature, and 
to our knowledge have not been found experimentally. The lowest value of L reported in 0] is 0.43, which is even 
larger than the value of Fig. |]. However, 0.0625 would easily be reachable by increasing the jet radius. 
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